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Abstract 

The free energy density, or pressure, of QCD has at high temperatures an expansion in the 
coupling constant g , known so far up to order g 5 . We compute here the last contribution which 
can be determined perturbatively, g®ln(l/g), by summing together results for the 4-loop vacuum 
energy densities of two different three-dimensional effective field theories. We also demonstrate that 
the inclusion of the new perturbative g 6 ln(l /g) terms, once they are summed together with the so 
far unknown perturbative and non-perturbative g 6 terms, could potentially extend the applicability 
of the coupling constant series down to surprisingly low temperatures. 

PACS numbers: ll.10.Wx, 05.70.Ce, ll.15.Bt, 12.38.Bx, 12.38.Mh 


1 


I. INTRODUCTION 


Due to asymptotic freedom, the properties of QCD might be expected to be perturbatively 
computable in various “extreme” limits, such as high virtuality, high baryon density, or high 
temperature. We concentrate here on the last of these circumstances, that is temperatures T 
larger than a few hundred MeV. 

The physics observable we consider is the pressure, or minus the free energy density, 
of the QCD plasma. Potential phenomenological applications include the expansion rate 
of the Early Universe after it has settled into the Standard Model vacuum, as well as the 
properties of the apparently ideal hydrodynamic expansion observed in on-going heavy ion 
collision experiments, just shortly after the impact. 

In these environments, it turns out that the naive expectation concerning the validity of 
perturbation theory is too optimistic. Indeed, even assuming an arbitrarily weak coupling 
constant g , perturbation theory can only be worked out to a finite order in it, before the 


0 . 


serious infrared problems of finite temperature held theory deny further analytic progress 
For the pressure, the problem is met at the 4-loop order, or 0(g e ). 

This leads to the interesting situation that there is a definite limit to how far perturbation 


0 . 


theory needs to be pushed. So far, there are known loop contributions at orders 0(g 
0(9 3 ) Q, 0(9-111(1/9)) 0 , 0(9-) 0, and 0(g 5 ) 0 . There is also an all-orders numeri¬ 
cal result available for a theory with an asymptotically large number of fermion flavors j^|. 
The purpose of the present paper is to collect together results from two accompanying pa¬ 
pers |9| |10j] , allowing to determine analytically the last remaining perturbative contribution, 
0(g 6 ln(l/g)), for the physical QCD. 

It must be understood that even if computed up to such a high order, the perturbative 
expansion could well conver ge o nly very slowly, requiring perhaps something like T TeV, 


to make any sense at all 


/erge only 

mu 


With one further coefficient available, we can to some 
extent now reinspect this issue. To do so we actually also need to assume something about 
the unknown 0(g e ) term, since the numerical factor inside the logarithm in (D(g 6 ln(l/c/)) 
remains otherwise undetermined. Therefore, our conclusions on this point remain on a 
conjectural level, but turn out to show nevertheless a somewhat interesting pattern, which 
is why we would like to include them in this presentation. 

Finally, it should be stressed that even if the perturbative expansion as such were to 
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remain numerically useless at realistic temperatures, these multiloop computations are still 
worthwhile: the infrared problems of finite temperature QCD can be isolated to a three- 
dimensional (3d) effective held theory and studied non-perturbatively there with simple 
lattice simulations However, to convert the results from 3d lattice regularisation to 3d 
continuum regularisation, and from the 3d continuum theory to the original four-dimensional 
(4d) physical theory, still necessitates a number of perturbative “matching” computations. 
Both of these steps are very closely related to what we do here, although we discuss explicitly 
only the latter one. 


II. THE BASIC SETTING 


We start by reviewing briefly how it is believed that the properties of QCD at a finite 
temperature T can be reduced to a number of perturbatively computable matching coeffi¬ 


cients, as well as some remaining contributions from a series of effective held theories 


Our presentation follows mostly that in 


Q, 


Q. 


but there are a few significant differences. 

The underlying theory is finite temperature QCD with the gauge group SU(iV c ), and 
Nf havors of massless quarks. In dimensional regularisation the bare Euclidean Lagrangian 
reads, before gauge fixing, 

rPh 

1 0 


5, 


QCD 


dr /d d xL 


QCD, 


C 


= T KvKv + ^luDnlp , 


QCD - ^ M fllS 1 Y I 


( 2 , 1 ) 

( 2 . 2 ) 


where (5 — T d = 3-2e, ft, v = 0 = d^A“- d„A°+ gf° b °A b li A%, £> p = cl, 

A^ — A°T n , "(2 — 7 ,,, { 7^,7 u} — 25^, and carries Dirac, color, and flavor indices. 

Denoting the generators of the adjoint representation by ( F a )b c = —i f abc , we define the 
usual group theory factors, 


CA^ab 

= [F c F c ] ab , 

C F 5ij 

= [T a T a ] ij: 

(2.3) 

T A S ab 

= Tr F a F b , 

T F 5 ab 

= Tr T a T b , 

(2.4) 

d A 

= S aa = N 2 C - 1 , 

d F 

= fiii = T F dA/C F . 

(2.5) 


Obviously Th = Ca■ For the standard normalisation, with Nf quark havors, Ca = N c , Cp = 
(iV c 2 - 1)/(2iV c ), T a = N c , T f = N f / 2, d A = N/ - 1, d F = N c N f . 
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We use dimensional regularisation throughout this paper. The spatial part of each mo¬ 
mentum integration measure is written as 


f = [ _ —2e 

l p -J (27 T) d A 


p 2e 


/ e 7 \ e f d d p 


\ 47 tJ J (2n) d . 


( 2 . 6 ) 


where /i = /i(e 7 /47r) 1 / 2 , and the expression in square brackets has integer dimensionality. 
From now on we always assume implicitly that the factor p~ 2e is attached to some rele¬ 
vant coupling constant, so that the 4d g 2 is dimensionless, while the dimensionalities of 
g E , A® and g^ , to be introduced presently, are GeV. 

The basic quantity of interest to us here is minus the free energy density /qcd(^), or the 
pressure Pqcd(T), defined by 


Pqcd (T) = (7 In j VA“ P^CV>exp(-is Q cD), 


(2.7) 


where V denotes the d -dimensional volume. Boundary conditions over the compact time¬ 
like direction are periodic for bosons and anti-periodic for fermions. Moreover, we assume 
Pqcd (T) renormalised such that it vanishes at T = 0. To simplify the notation, we do not 
show the infinite volume limit explicitly in the following. 

At high temperatures and a small coupling, there are parametrically three different mass 
scales in the problem, ~ 2tcT, gT , g 2 T |l3| . All the effects of 
be accounted for by a method called dimensional reduction 


lie hard mass scale ~ Znl can 
Specifically, 


jllc Hell 1 

3Q. 


Pqcd(T) = p E (T) + | hi J VAl VA^ exp(-5' E ), (2.8) 

<Se = Jd d xC E , (2.9) 

C E = ^TrFZ l + Tr[D k ,A 0 ] 2 + m 2 E TrAl + \ ( £\TrA 2 0 ) 2 + \ ( g ) TrA* + ... . (2.10) 


Here k — 1, F ^ = (i/g E )[Dk, Di], Dk = dk — ig E Ak, and we have used the shorthand 
notation = A^T a , A 0 = Agf a , where T a are Hermitean generators of SU(1V C ) normalised 
such that Tr T a T b = S ab /2. Note that the quartic couplings Ag\ A® are linearly independent 
only for N c > 4. 

The relation in Eq. (|ZHD contains hve different matching coefficients, p E , m El g E , Ag\ Ag\ 
We are interested in the expression for Pqcd (7") up to order 0(g e T 4 ). They will then have 
to be determined to some sufficient depths, as we will specify later on. Let us here note that 
the leading order magnitudes are p E ~ T 4 , m| ~ g 2 T 2 , g\ ~ g 2 T, A^ ~ g 4 T, A^ ~ g A T. 
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Apart from the operators shown explicitly in Eq. eh, there are of course also higher 


order ones in C E . The lowest such operators have been classified in Their general 

structure is that one must add at least two powers of D k or gA 0 , to the basic structures 
in Eq. ehh Since higher order operators are generated through interactions with the scales 
that have been integrated out, ~ 2nT, they must also contain an explicit factor of at least 
g 2 . For dimensional reasons, the schematic structure is thus 


SC, 


DuDi 


rA 


( 2 . 11 ) 


3 (2ir Tf 

To estimate the largest possible contributions such operators could give, let us assume the 
most conservative possibility that the only dynamical scale in the effective theory is r-u gT. 
By dimensional analysis, we then obtain a contribution 

Spqcd(T) _ , r 2 (gT ) 2 


SCt- 


-{gTf ~ g 7 T 3 


( 2 . 12 ) 


T ~~~ " (2vrT) 2 ' 

Therefore, all higher dimensional operators can be omitted from the action in Eq. EH , if 
we are only interested in computing pqcd(T) up to order (T(g 6 T 4 ). 

The theory in Eq. EH contains still two dynamical scales, gT , g 2 T. All the effects of 

Specifically, 

(2.13) 

(2.14) 

(2.15) 


the “color-electric” scale, pT, can be accounted for by integrating out Aq 

^ hi J VA a k VA a 0 exv(-Sv) = p u (T) + ^ In j VA a k exp(-A M ), 


Am = d d xC 


M, 


£m — - Tr + ... , 


where F ki = ( i/g M )[D k , A], D k = d k - ig M A k , and A k = AlT a . 

The relation in Eq. EH) contains two matching coefficients, Pm^m, which again have 
to be determined to sufficient depths. At leading order, p M ~ m|T, g ^ ~ g\. In addition, 
there are also higher order operators in Eq. EH) - The lowest ones can be obtained by 
imagining again that we apply at least two covariant derivatives to Eq. EH , together with 
at least one factor g\ brought in by the interactions with the massive modes. This leads to 
an operator 

5C M ~ g 2 J^C M . (2.16) 

m E 

The only dynamical scale in the effective theory being ~ g 2 T, dimensional analysis indicates 
that we then obtain a contribution of the order 

^PQCd(T) 2 (g 2 T) 2 {„2rp^3 ^9rp3 ( 0 i7( 

———~ (h-M ~ 9 E — 3 —yg 1 ) ~ g 1 ■ ( 2 - 17 ) 

l m E 
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Therefore, higher dimensional operators can again be omitted, if we are only interested in 
the order 0(g 6 T 4 ) for pq C d(T). 

After the two reduction steps, there still remains a contribution from the scale g 2 T, 

Pg(T) = 2T4£exp(—S M ) , (2.18) 

with S'm in Eqs. firm cnsD. Since £m only has one parameter, and it is dimensionful, 
the contribution is of the form 

Pg(T) ~ Tg 6 M . (2.19) 

The coefficient of this contribution is, however, non-perturbative [l[ Q] . 

In the following sections, we proceed in the opposite direction with regard to the presen¬ 
tation above, from the “bottom” scale g 2 T , producing p G (T), through the “middle” scale 
gT , producing pu(T), back to the “top” scale 27rT, producing pe(T). We collect on the way 
all contributions up to order g 6 T 4 , to obtain Pqcd(T) = Pe(T) + Pm{T) + p G (T). 


III. CONTRIBUTIONS FROM THE SCALE g 2 T 


The contribution to Pqcd(^) from the scale p ~ g 2 T is obtained by using the theory C M 
in Eq. (12.151) in order to compute p G (T), as defined by Eq. (12.181) . 

As is well known [1, Q, the computation involves infrared divergent integrals, starting 
at the 4-loop level. This is a reflection of the fact that Cm defines a confining held theory. 
Therefore, p G (T) cannot be evaluated in perturbation theory. 

What can be evaluated, however, is the logarithmic ultraviolet divergence contained in 
p G {T). For dimensional reasons, the non-perturbative answer would have to be of the form 


Pg(T) i r3 ,9m 


Tp 


—2e 


(47t) 4 . 


a G (- + 8ln-^-)+f3 G + 0(e) 
2mM y 


(3.1) 


where tum = CaP m- Now, because of the super-renormalisability of Cm, the coefficient a G 
can be computed in 4-loop perturbation theory, even if the constant part (3 G cannot B- 
Of course, if we just carry out the 4-loop computation in strict dimensional regularisation, 
then the result vanishes, because there are no perturbative mass scales in the problem. This 
means that ultraviolet and infrared divergences (erroneously) cancel against each other. 
Therefore, we have to be more careful in order to determine a G . 
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To regulate the infrared divergences, we introduce by hand a mass scale, m^., into the 
gauge held (and ghost) propagators. This computation is described in detail in |9j. Individ¬ 
ual diagrams contain then higher order poles, like 1 /e 2 , as well as a polynomial of degree up 
to nine in the gauge parameter £. However, terms of both of these types cancel in the final 
result, which serves as a nice check of the procedure. 

As a result, we obtain 

Pg(T) 


T/r 


—2e 


j r , ‘ i 

dACA l 4 ^ 


«g(- + 8 In 7 ^-) +^ G (0 + C> ( e ) 


(3.2) 


'6 2?7Tg / 

where is used to denote that only the coefficient a G multiplying 1 /e is physically 

meaningful, as it contains the desired gauge independent ultraviolet divergence, defined 
in Eq. (EU>. The value of the coefficient, obtained by extensive use of techniques of sym¬ 
bolic computation (implemented []jJ in FORM 

43 157 


LiCU. uy 

Q), is D 


“ G= % 


6144 


-7T 2 


0.195715 . 


(3.3) 


On the contrary, the constant part $g{£) depends on the gauge parameter £, because the 
introduction of m G breaks gauge invariance, and has nothing to do with /3q in Eq. (EH). 

IV. CONTRIBUTIONS FROM THE SCALE gT 


We next proceed to include the contribution from the scale gT, contained in Pm CO, as 
defined by Eq. (12.131) . 

By construction, Eq. EUD assumes that all the infrared divergences of the expression on 
the left-hand-side are contained in vn(T ). dehned in Eq. (12.181). and determined in Ea. EH- 
Therefore, if we compute the functional integral (T/V) ln[/ VA^VAq exp(— S'e)] using strict 
dimensional regularisation (i.e., without introducing by hand any mass ?»g for the gauge held 
Aj), whereby Pg(T) vanishes due to the cancellation between infrared and ultraviolet diver¬ 
gences mentioned above, we are guaranteed to obtain just the infrared insensitive matching 


coefficient Pm(0- This is exactly the computation we need, and carry out in 10, [19]. It 


may be mentioned that we have checked explicitly the infrared insensitivity of the result, by 
giving an equal mass to both A 0 and Aj in the 4-loop expression for the functional integral, 
and then subtracting the graphs responsible for p G (T), with the same infrared regularisation. 
This result is also independent of the gauge parameter. 
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Keeping terms up to order 0(g e T 4 ), the full outcome for pm(T) is 


Pm(T) 

2 >- 2e 


(An 

1 

(47t) 2 

1 

(47t) 3 

1 

(47t) 4 

1 


■ d A m1 - + 0(e) 


d A CAg E m E 


d A C A g E m E 


1 3 

4e _ 4 
1 


In 


/i 


2 m\ 


+ 0(e) 


11 


—--tt 2 + —ln2 + 0(e) 
24 6 6 W 


d A C\g E 


«M 


(—8 In —) + ( 3 m + @( e ) 

M 2mE y 


(An) 


d A (d A + 2)X ( E ) ml 


where 


0 


(An 


1 A 2 ^ — 1 x ( 2)_2 
mUA Ar TYIt? 


N r 


-4+°W 

-i+0( £) 


43 4Q1 

«m = — — . 7T 2 Rj 0.555017 . 


(4.1) 


(4.2) 


32 6144 

The finite constant /3m can be expressed in terms of a number of finite coefficients related 
to 4-loop vacuum scalar integrals Jjioj], but we do not need it here. 

In addition to pm(T), we also need to specify the effective parameter g ^ appearing in 
Cm, to complete contributions from the scale gT. It is of the form 


9 m ^ 9 e { 1 + 0 (gl/m E )) , 

where the next-to-leading order correction is known (see, e.g., 0), 


(4.3) 

but not needed here. 


V. CONTRIBUTIONS FROM THE SCALE 2nT 


The contributions from the scale 2nT are contained in the expressions for the parameters 
of the previous effective theories, as well as in p E (T). We write these as 


/i 2t p E (T) — T 4 q?ei + g 2 (oiE2 + 0(e)^j 


2 ATI 2 

m E = T 


9 e = 


T g 2 + («E7 + /3 E3 6 + C^e 2 )) + 0(g 6 ) 


Ag> = T 


-(47t) ; 


■(/ 3 E4 + 0(e))+0(g 6 


(5.1) 


+ jhf (" es + + jhf (^ ei + ° + 

fi ,2 («E4 + CtE5 e + C^(C)) + ( a E6 +/3 e 2C + C^(c 2 )) + ^(l/ 6 ) j (5-2) 


(5.3) 

(5.4) 
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A<?> = T 


+ 0(e))+0(g 6 


(5.5) 


where g 2 is the renormalised coupling. We have named explicitly («e, Pe) the coefficients 
needed up to order 0(g 6 ). The actual values for those needed at order 0(g b ln(l/c/)), denoted 
by ctE, are given in Appendix El The additional coefficients needed at the full order 0(g 6 ) 
are denoted by Pe] some of these are also known (for Pea , Pea , e.g., see | 3 ])- The rest of the 
terms contribute only beyond 0(g b ). 

The expression for Pe(T) is simply the functional integral in Eq. (12.71) . calculated to 4-loop 
level in the MS scheme, but without any resummations. The only physical scale entering 
is thus 2ttT. The calculation has so far been carried out only to three loops jfl [ll| so that 
Pei is not known. Even when performed with the fully renormalised theory, the results in 
general contain uncancelled 1/e poles, as explicitly seen in the 3-loop expression in Eq. (I A. 31) 
for «E 3 - These only cancel when a physical fully resummed quantity is evaluated, i.e., in the 
sum Pqcd = Pe + Pm +Pg■ Similarly, m|, g|, Aj^ can be obtained for instance from suitable 
2-, 3-, and 4-point functions, respectively. 


VI. THE COMPLETE RESULT 


Combining now the results of Secs. urn Gum and expanding in g , we arrive at 
Pqcb(T) Pe(T) + pm{T) + pc(T) 


T 4 /i 


— 2e 


T 4 /i 


— 2e 


= 9° i«El 


+ 9 |«E2 

9 3 / d A 3/2 


(47t) l 3 

g A 


■a 


E4 


(4tt 
9 5 


2 i aE3 — d A C A 


13 g \ 1 

a E 4 1 --1“ T + l n “—~ f/2 ) + 7 OL E 5 

TT^ A ' 


(47t) 3 
^6 


d A Ol E4 


4e 4 2 gTo& 

1 2 789 7T 2 11 0 

— c^E 6 — C7 a ( — H-— — In 2 

2 6 A V24 6 6 


J^y{(dEi ~^d A a E4 


( d A + 2)/3e4 H- jy — Pe5 


-d A C A 


— («E6 + «E5«E7 + 3«E4«E7 + Pe2 + «E4/3es) 




+ («E6 + aE4«E7) ( --h In i/o ! 

V y V4e 2 //T«e / 4 /J 
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+(1aC\ 

■0(g 7 ) + 0(e 


P m + Pg + a M ( —h 8 In- 1/2 

e 2 gTc^l 


+ « G ( - + 


; In 


A_V 


2g 2 TC A J\ 


( 6 . 1 ) 


Utilising the expressions in Appendix El the terms up to order 0{g 5 ) reproduce the known 
result in jj]. 

For the contribution at order 0(g 4 ), the 1/e divergence in a E3 (cf. Eq. (IA.3I) 1 and 
the 1/e divergence from pm(T), shown explicitly in Eq. (16.1(1 . cancel. This must happen 
since Pqgd(T) is a physical quantity. The associated JV s also cancel, but a physical effect 
ln[mE/(27rT)] ~ \n{ga^) remains 5], 

For the contribution at order 0(g 6 ), a number of unknown coefficients remain (the /?e’s, 
(3m, (3g), but a similar cancellation is guaranteed to take place. In addition, the result must 
be scale independent to the order it has been computed. The first point can be achieved by 
Pei (the other /?e’s are finite), so that it has to have the structure 


Pei = gUCU(cKE6 + OlEAOlEl)— ~ (IaC\{oi m + Oto)- + pE6 , 


( 6 . 2 ) 


where Peq does not contain any 1/e poles. The latter point can be achieved by adding and 
subtracting ln[/Z/(27rT)]’s, such that J1 gets effectively replaced by 2tcT in the logarithms 
visible in the 0(g 6 ) term in Eq. (16.11) . The ln[/Z/(27rT)]’s left over, together with those 
coming from the /3 e ’s, serve to cancel the effects from the 2-loop running of g 2 (fl) and 1-loop 
running of g 4 (fi) in the lower order contributions, without introducing large logarithms. 

This general information is enough to fix the contributions of order 0(g 6 ln(l/g)) to 
Pqcd(T). Indeed, after inserting Eq. (16.21) and reorganising the logarithms appearing in the 
Pe’s as mentioned, there remains a logarithmic 4-loop term, 


Pqcd(T) 

T 4 /i- 2e 


g e ln(l /g) 


= 9 


: 3jaCa /, 
(47r) 4 


( («E6 + aE4<*E7) lii(jaE4 2 ) “ o M + 2a G In (gC 1 / 2 ^ J 


(6.3) 


where «e 4 is in Ecp (TOl) . cx E 6 is hi Eq. (IA.6I) . cx E 7 is in Eq. (ED, ckm is in Eq. (ED, and 
ckg is in Eq. (ED- Note that there are logarithms of two types, with different non-analytic 
dependences on group theory factors inside them. Eq. ED is our main result. 


Following 



, let us finally insert N c 


3, and give also the numerical values for the 
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various coefficients, for an arbitrary Nf. We obtain 


where 


Pqcd(T) = 


8n 2 

~45 


rj~i4 


ST' f a s{p )\ 1 / 2 

S *><(—) 


i=0 


Po 

Pi 

P-2 

P3 

P4 


21 Ar 

1 H- N f , 

32 1 ’ 


= 0, 


30(1 + 1 N ,f\ 

1 QC 1 1 

237.2 + 15.96IV/ - 0.4150 N 2 + — (l + - N f ) ln[— (l + - N f 


165 


( 1 + 4jv / )( i -|jv / ) i „ 


12 

1/2 


P_ 
2t tT ’ 


p 5 = (l + - N f ) -799.1 - 21.967V/ - 1.926 N 2 f 


495 


0 + 6^) 0-55^) 111 


P 


33 


2ttT 


P6 = 


-659.2 - 65.89 N f - 7.653iV; 

1 

6 




P 


2t rT 


In 


EK*) 


-475.6 In ^ + ?,(«/) In 2 ^ + q b (N f ) 111 ^ + qt(JV/) 


(6.4) 


(6.5) 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 


( 6 . 11 ) 


where q a (Nf ), %(Nf ), q c (Nf) are a s -independent polynomials in Nf. Two of them, q a (Nf), 
qb(Nf), can already be written down because they just cancel the //-dependence arising from 
the terms of orders a s (fi),a 2 (fi ): 

w 2 \ 2 

Vi-11 I- Nf 

12 

(6.13) 


^) = -“(i4^)(i-|<. 


16 l 12 33 

%(Nf) = 2932.9 + 42.83/V/ - 16.48/Vj + 0.2767V) 


The third one, q c {Nf), remains however unknown. 


VII. THE NUMERICAL CONVERGENCE 

This Section is devoted to a numerical discussion of the result. Since the 0(g 6 ln(l/g)) 
term cannot be given an unambiguous numerical meaning until the 0(g 6 ) term is specified, 
we have to present the result for various choices of the latter. In the relevant range of T/Ayjg 


11 



the outcome will depend sensitively, even qualitatively, on this uncomputed term. One choice 


will be seen to agree with 4d lattice data down to about T/A 


MS 


2...3. Since however 


dimensional reduction, that is an effective description of QCD via the theory in Eq. (inn , 
is known to break down at about this point, and we have only kept a finite number of terms 
in the expansion following from Eq. mm), this cannot really be considered a prediction, 
even if the eventual computation of the 0(g 6 ) term gave just the appropriate value. It is 
just an observation that a smooth transition from the domain of validity of our results to a 
domain of different approximations should be possible. 

A standard procedure in the discussion of perturbative results would be to take the 
expansion in Eq. (ED and to study whether its scale dependence is reduced when further 
orders of perturbation theory are included. As is well known since jfij], this fails for the 
pressure, unless T A^g. Related to this, the numerical convergence of the perturbative 


expansion is known to be qui 
below the electroweak scale 17 



oor for any fixed scale choice, at least for temperatures 


11 , Va 


The new term we have computed does not change 
this general pattern. But the culprit is known: it is Pm(T) + pc(T) emerging from the 3d 
sector of the theory, where the expansion parameter is only g^/(nmE) ~ g/i r. In contrast, 
for Pe(T) as well as for, say, jet physics, the expansion parameter is a s /n, and there are 
good reasons to expect numerical convergence to be much better. 

For these reasons, we will only discuss the sensitivity of the result on the so far unknown 
0(g 6 ) coefficient, as well as the slow convergence of the 3d sector, in the following. For 
simplicity, we only consider the case N c = 3, Nf = 0 here. 

As in |l4j, the actual form we choose for plotting contains Pm{T) + Pq(T) (Eqs. (14.111 + 
(13.11) 1 in an “un-expanded” form, that is, with the , ,q| inserted from Eqs. f|5,2H . (15.3(1 . and g ^ 
from Eq. mi). This means that we are effectively summing up higher orders: the C>(g 3 )-term 
is really 0(g 2 + p 4 ) 3//2 , while the 0(g 6 ln(l/g)) term contains a resummed coefficient, being 
then effectively 0((g 2 + q 4 ) 3 lnfl/olb We proceed in this way because then a comparison 
with numerical determinations 14 { of the slowly convergent part Pm(T) + Pg(T) is more 
straightforward, and also because the resummations carried out reduce the /7-dependence of 
the outcome. However, we have checked that the practical conclusions remain the same even 
if we plot directly the expression in Eqs. (I(1.4l) - (l(i.lll) (but with a larger scale dependence). 

To be specific, the genuine 0(g 6 ln(l /g) + g 6 ) contribution, which collects the effects from 
all the terms involving the /3 E ’s, /3m, As, «m, and etc hi Eq. m, is now written in the form 
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FIG. 1: Left: perturbative results at various orders (the precise meanings thereof are explained 
in Sec. m , including 0(g 6 ) for an optimal constant, normalised to the non-interacting Stefan- 
Boltzmann value psb- Right: the dependence of the 0(g e ) result on the (not yet computed) 
constant, which contains both perturbative and non-perturbative contributions. The 4d lattice 
results are from 


(specific for N c = 3, Nf = 0, where g\ r-u i /g), 


Pqcd(T) 

T/i 


—2 e 




(ctm T 2q:g) In —5—b <5 
9e 


(7.1) 


g 6 ln(l/ g) 

while the remaining 0(g 6 ) terms of Eq. (16.11) are contained in the resummed lower order 
contributions. The results are shown in Fig. ^ for various values of 6. The power of g 
labelling the curves indicates the leading magnitude of each resummed contribution. The 
scale is chosen as jl ~ 6.7T, as suggested by the next-to-leading order expression for g\ jl3 | . 
We observe that for a specific value of 5, the curve extrapolates well to 4d lattice data. 

While Fig. Q] looks tempting, the question still remains whether the good match to 4d 
lattice data with a specific value of the constant is simply a coincidence. This issue can be 
fully settled only once the constant is actually computed. However, we can already inspect 
how the slowly convergent part of the pressure, Pm+Pg, really behaves. 

The different finite terms in (pm + Pg) / {T g^) are plotted in FigEl The A {^-contributions 
are negligible. The results depend then essentially only on m|/ g^, which for N c = 3, Nf = 0 
is ~ 0.32 log 10 (T/Ajyjg) + 0.29. We observe that the leading 1-loop term 0(g 3 ) is 
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FIG. 2: The absolute values of the various terms of the slowly convergent expansion for Pu(T) + 
Pg{T), normalised by Tg^. 

dominant for T/A^g > 10, the 3-loop term 0(g 5 ) is rather big, bigger in absolute value than 
the 2-loop term 0(g A ) within the T-range of the figure, while the 4-loop term is always very 
small. Therefore, while it is well possible that there is again a big “odd” 0(g 7 ) contribution, 
it is perhaps not completely outrageous either to hope that the convergence could also 
already be reasonable, once the full 0(g 6 ) contribution is included. If this were the case, 
then all higher order contributions would have to sum up to a small number. 

Finally, it is perhaps interesting to remark that at the time of the numerical lat¬ 
tice Monte Carlo study in ref. [ffij], nothing was known about the coefficient /3ei, which 
was therefore set to zero (cf. Eq. (4) in [14;]), while the part Pm(T) + Pq{T) was deter¬ 
mined non-perturbatively. But this means that a logarithmic term coming from the scale 
27 tT, ~ —g 6 (a m + a g) hi[/i/(27rT)], was missed. With the scale choice p = p E = g\ 
within results obtained with this converted to a missing 0{g & \n{l/g)) contribution 
g 6 (2a M + 2 <ug) ln(l /g). With the same scale choice the non-perturbative part, on the other 
hand, contributed — g & otu hi(l/ < 7 ) and led to the wrong curvature of the pressure seen at 
small T/A^jg. Adding the missing part, which now has been computed, leads to a total of 
g 6 (a m + 2<ug) ln(l /g) with the opposite sign and the correct (he., the one seen in 4d lattice 
measurements) curvature in Fig. Q (for small values of <5). Therefore the 0(yg & \n(l/g)) terms 
are indeed physically very relevant. 
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VIII. CONCLUSIONS 


We have addressed in this paper the 4-loop logarithmic contributions to the pressure of 
hot QCD. Physical (regularisation independent) logarithms can only arise from a ratio of two 
scales. Since there are three parametrically different scales in the system, 2ttT, gT, g 2 T, there 
are then various types of perturbatively computable logarithms in the 4-loop expression for 
the pressure: 


1. Logarithms of the type g 6 ln[(27rT)/(g 2 T)]. The coefficient of these is computed in Q, 
and given in Eq. m- 


2 . 


Logarithms of the type g e ln[(2irT)/(gT)\. The coefficient of these is computed in 
and given in Eq. 



3. Logarithms related to the running of the coupling constant in the 3-loop expression 
of order 0(g A h\[(2nT) / (gT)]). Their MS coefficient can be seen in the first term 
in Eq. o, but it depends on the scheme, and can in principle even be chosen to 
vanish. 


Logarithms of the first and second types can be written in many ways: it may be more 
intuitive, for instance, to reorganise them as 

ff°G ln (“^) + 9 6a M ln (“^") = 9 e ( aM + "g) hl (^“) + 9&aG hl (^r) ' ( 8 ' 1 ) 

The existence of three kinds of logarithms is somewhat specific to non-Abelian gauge 
theory. In QED, in particular, none of the logarithms appear. This is due to the fact 
that the effective theories we have used for their computation, Eqs. HUB, (I215D , are non¬ 
interacting (apart from a term K in Eq. firm which does not lead to logarithms). 
Therefore we have nothing to add to the known 0(g 5 ) QED result obtained in [23|. In the 
0 4 scalar theory, on the other hand, there is a logarithm of the second type, and also one 
somewhat analogous to the third type. Their coefficients were already computed in Q. 

There are interesting checks that can be made on the various logarithms mentioned, using 


methods completely different from those employed here. For instance, logarithms of t 
and second types could in principle be seen with 3d lattice Monte Carlo methods [ikl 
well as with stochastic perturbation theory 


le first 
2 fi| . as 


A very interesting analytical check would 
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be to compute the 4-loop free energy directly in 4d in strict dimensional regularisation, 
but without any resummation. By definition, this computation produces the coefficient 3e\ 
in Eq. m [ll| , and one check is that the result must contain the 1/e divergences shown 
in Eq. (lb.21) . 

To complete the free energy from the current level 0(g 6 ln(l/g)) to the full level 0(g 6 ), 
would require significantly more work than the computation presented here. More specif¬ 
ically, there are contributions from all the scales in the problem, ranging from 2irT (the 
coefficients /?ei, Pes), through gT (the coefficient /3m), down to the non-perturbative 
scale g 2 T (the coefficient /3q). This then requires carrying out 4-loop finite temperature 
sum-integrals, 4-loop vacuum integrals in d = 3 — 2e, 4-loop vacuum integrals in 3d lattice 
regularisation, and lattice simulations of the pure 3d gauge theory in Eq. (I?~T51) . 

Nevertheless, given the potentially important combined effect of all these contributions, 
as indicated by Fig.|T| such computations would clearly be well motivated. 
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APPENDIX A: MATCHING COEFFICIENTS 


In Eqs. dn}-(ESD we have defined a number of matching coefficients, the cte’s and /?e’s. 


For the oe’s, the following expressions can be extracted from HQQ: 
_2 

■(4c l a + 7 dp'j , 
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+ T ln2 


(A.3) 
(A.4) 
(A.5) 

(A.6) 
(A.7) 


Note that with our notation, the 1-loop running of the renormalised coupling constant goes 


as 


/O) = 9 2 {V 0 ) 


^ (HC'a — 4T f ) 


^ 4 (/^o) 

(47r) 2 
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